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Abstract. We study the asymptotic behavior of the wave func- 
tion in a simple one dimensional model of ionization by pulses, 
in which the time-dependent potential is of the form V(x, i) = 
—2S(x)(l — e~~ xt cosui), where S is the Dirac distribution. 

We find the ionization probability in the limit t — > oo for all A 
and u>. The long pulse limit is very singular, and, for u> = 0, the 
survival probability is const A 1 / 3 , much larger than O(A), the one 
in the abrupt transition counterpart, V{x, t) — 6(x)l{ t >i/\y where 
1 is the Heaviside function. 



1. Introduction 

Quantum systems subjected to external time-periodic fields which 
are not small have been studied in various settings. 

In constant amplitude small enough oscillating fields, perturbation 
theory typically applies and ionization is generic (the probability of 
finding the particle in any bounded region vanishes as time becomes 
large). 

For larger time-periodic fields, a number of rigorous results have 
been recently obtained, see [5] and references therein, showing generic 
ionization. However, outside perturbation theory, the systems show a 
very complex, and often nonintuitive behavior. The ionization fraction 
at a given time is not always monotonic with the field [Bj. There 
even exist exceptional potentials of the form S(x)(l + aF(t)) with F 
periodic and of zero average, for which ionization occurs for all small 
a, while at larger fields the particle becomes confined once again |lUj . 
Furthermore, if 5(x) is replaced with smooth potentials f n such that 
f n —>■ 5 in distributions, then ionization occurs for all a if n is kept 
fixed. 

Numerical approaches are very delicate since one deals with the 
Schrodinger equation in W 1 x M + , as t — > oo and artefacts such as 
reflections from the walls of a large box approximating the infinite do- 
main are not easily suppressed. The mathematical study of systems in 
various limits is delicate and important. 
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In physical experiments one deals with forcing of finite effective du- 
ration, often with exponential damping. This is the setting we study 
in the present paper, in a simple model, a delta function in one dimen- 
sion, interacting with a damped time-harmonic external forcing. The 
equation is 

(1) i ?t = ( - ^ - 25(x) (1 - A(t) cosM)) ) 

where A(t) is the amplitude of the oscillation; we take 

(2) = ^(0, x) G Cq 00 ; A(t) = ae~ xt ; a = 1 

(The analysis for other values of a is very similar.) The quantity of 
interest is the large t behavior of ip, and in particular the survival 
probability 

(3) P B = lim P(t, B) = lim / \ip(t,x)\ 2 dx 

t^oo t— >oo J g 

where B is a bounded subset of R. 

Perturbation theory, Fermi Golden Rule. If a is small enough, P 
decreases exponentially on an intermediate time scale, long enough so 
that by the time the behavior is not exponential anymore, the survival 
probability is too low to be of physical interest. For all practical pur- 
poses, if a is small enough, the decay is exponential, following the Fermi 
Golden Rule, the derivation of which can be found in most quantum 
mechanics textbooks; the quantities of interest can be obtained by per- 
turbation expansions in a. This setting is well understood; we mainly 
focus on the case where a is not too small, a toy-model of an atom 
interacting with a field comparable to the binding potential. 

No damping. The case A = is well understood for the model ([1]) in 
all ranges of a, see [13]. In that case, P(t, A) ~ t~ 3 as t —>■ oo. 

However, since the limit A —>■ is singular, little information can be 
drawn from the A = case. 

For instance, if u = 0, the limiting value of P is of order A 1//3 , while 
with an abrupt cutoff, A(t) = l{t:t<i/A}) the limiting P is 0(A) (as 
usual, Is is the characteristic function of the set S). 

Thus, at least for fields which are not very small, the shape of the 
pulse cut-off is important. Even the simple system ([T]) exhibits a highly 
complex behavior. 

We obtain a rapidly convergent expansion of the wave function and 
the ionization probability for any frequency and amplitude; this can be 
conveniently used to calculate the wave function with rigorous bounds 
on errors, when the exponential decay rate is not extremely large or 
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small, and the amplitude is not very large. For some relevant values of 
the parameters we plot the ionization fraction as a function of time. 

We also show that for uo = the equation is solvable in closed form, 
one of the few nontrivial integrable examples of the time-dependent 
Schrodinger equation. 

2. Main results 

Theorem 1. Let if)(t,x) be the solution of (TJP with initial condition 
ij G L7 °°. Let 

(1) 9m,n = 9m,n(°) = l f e~^^^(x)dx' 

Then as t — ► oo we have 

(2) ij)(t, x) = r(A, u)e u e- lxl (l + t~ 1/2 h(t, x)) 
where \h(t, x)\ < C , Va; G M., Vt G M + ; and where 

(3) r(A, u) = [-Ax,-! - A hl + 2(/ ,oUi 



where A m ^ n = A m ^ n (a) solves 

(4) (Va + ncu - imX - i)A m>n = -~A m+ltn+1 - -A m+1 , n -x + g mjTi 
There is a unique solution of Oj satisfying 



(5) ^(l + W)ie- 6 V^H|A mi „|<oo 

m,n 

where b > 1 is a constant. It is this solution that enters (T3|). 

There is a rapidly convergent representation of r(A, u), see §3.51 
Clearly, |r(A,a;)| 2 is the probability of survival, the projection onto 
the limiting bound state. 

2.1. u = 0. 

Theorem 2. (i) For u = we have 
(6) 

r ( A )=/ 1 I , - n / <?(*0exp j=— A 2 



exp W e " P 2 ( -i + e -p)Vsr(pA)3/» dp ) dfaip 



" OO 
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where g(k) = g k ,o- 

(ii) We look at the case when ipo = e~' x ' ; the bound state of the 
limiting time-independent system. Assuming the series ofr(X) is Borel 
summable in A for arg A G [0, |] (summability follows from (0|), but the 
proof is cumbersome and we omit it), as A — > we /iawe 

(7) r(A) ~ 2- 2 / 3 (-3 ? ) 1 / 6 7r- 1 / 2 r(2/3)e^t A 1 / 6 

Note: The behavior (jTJ) is confirmed numerically with high accuracy, 
constants included, see §5.31 

We also discuss results in two limiting cases: the short pulse setting 
(see $6]) and the special case A = (see ^7J). 

3. Proofs and further results 

3.1. The associated Laplace space equation. Existence of a strongly 
continuous unitary propagator for (TTJ) (see [26J v. 2, Theorem X.71) im- 
plies that for ip G L 2 (M d ), the Laplace transform 

POO 

#,p):= / ^{;t)e^dt 
Jo 

exists for 9ft(p) > and the map p — > i>(-,p) is L 2 valued analytic in 
the right half plane 

p G M = {z : > 0} 
The Laplace transform of (TTJ is 
d 2 



(1) 



dx 
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+ ip ifj(x,p) = iipo — 25(x)^(x,p) 



+8(x) (^(x, p — ico + A) + ip(x, p + icu + A) j 
Let p = ia + mA + inu and 

(2) y m ,n(%, a) — ^(x, ia + mX + inui) 
where ia G {z : < $sz < u, < < A}. 

Remark 3. Since the p plane equation only links values of p differing 
by mA + inu, m, n G Z, it is useful to think of functions of p as vectors 
with components m and n, parameterized by a. 

Thus we rewrite (pQ) as 

/ <9 2 \ 

(3) I ^ ~~ a ~ nuJ + J 



#o - 25(x)y m ^ n + (y m +i,„+i + y m +i : 



n— 1, 
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When |n| + \m\ 7^ 0, the resolvent of the operator 

d 2 . x 

— 7—17 + a + noo — imX 
ox 1 

has the integral representation 

(4) U m ,nf) (x) := G(K m>n (x-x'))f(x')dx' 
with 

K m) „ = y—ip = V o~ + ncj — imX 
where the choice of branch is so that if p G H, then K, mn is in the fourth 
quadrant, and where the Green's function is given by 

(5) G{K m>n x) = ^e-^-W 

Remark 4. If f(x) G Cq°, tzsmg integration by parts we have, as 
p — > 00 

,„s c(x) /r 

(/) ~ — + - 
P \P, 

where we regard g as an operator with p as a parameter; see also Remark 
Furthermore, ^ implies c(x) G L 2 . 

Define the operator € by 

(6) (Cy) [25(x)y TOin - (y m+ i,n+i + y m+ i )Tl -i)] 
Then Eq. ([3]) can be written in the equivalent integral form 

(7) y = igipo + <£y 
where q is defined in (]!]). 

Remark 5. Because of the factor K^ n in $5$, we have, with the iden- 
tification in Remark^ 



as p — > 00, for any function <p(p). 

3.2. Further transformations, functional space. In this section 
we assume ipo G C£j°. As in Remark HJ we obtain 

(8) itt fo = sM. + ( 1 



p \p S / 2 



for some C\{x) G L 2 . 
Let 

(9) h 1 (p) = h 1 (x,p) = c 1 (x)C(l m (t)) 
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For large p we have 

(10) Mp) = £lM + (_I_ 

Remark 6. As a function of x , h\(p) is clearly in L? and 

C^ihxip)) = Cl (x)l m (t) 

thus for t > 1 we have 

£-\h 1 (p)) = 

Substituting 

(11) y = y 1 + h 1 
in (0) we have 

(12) ?/! = ig^o - h + £ (hi) + €yi 

Let |/o = iQipo — hi + £ (/ii), then Remark [5] implies that for large p 

( 13 ) Vo = o (Jp 

and by construction y G L 2 as a function of x. 

We analyze CE2D in the space J% = L 2 (Z 2 x R, \\ ■ \\ b ), b > 1, where 



(14) ||y|| 6 := ^(l + |n|)le- 6 V^H|| y , 



12 

to, n|| i, 2 



We denote by the transformed wave function corresponding to y±. 
Writing y instead of yi, we obtain from (IT21) . 

(15) y = y + <ty 

Lemma 7. £ is a compact operator on 3%, and analytic in y/—ip. 

Proof. Compactness is clear since € is a limit of bounded finite rank 
operators. Analyticity is manifest in the expression of £ (see (jlj) and 

©)• □ 

Proposition 8. Equation 173]) /ias a unique solution iff the associated 
homogeneous equation 

(16) y = €y 

has no nontrivial solution. In the latter case, the solution is analytic 



in yja. 

Proof. This follows from Lemma [7] and the Fredholm alternative. □ 
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When m = 0, n = 0, and a = 0, £ is singular, but the solution is 
not. Indeed, by adding A > 0, to both sides of (j3J) we get the 

equivalent equation 

/ <9 2 \ 
( 17 ) ( g-j - ^ - w + imA + 1 [_a,a] I j/ m ,„ 

= iip + (1[-a,A] - 28{x)) y m , n + 8(x) (y m+ i, n +i + y m +x,n-x) 

Arguments similar to those when is absent show that the op- 

erator € associated to (fT7|) is analytic in y/a, thus y m n is analytic in 

3.3. Equation for A. Componentwise ([7j) reads 

(is) y y 

+ e" Km -" |x| [2|/ min (0) - (y m+ i,„ + i(0) +y m+1 ,„_i(0))] 



2k 



rn.n 



With A mjTl = j/ m , n (0), we have 

(19) (vct + - otiA - l)A m>n = --An+i.n+i - -A m+hn _x + g min 
where g m . n is defined in (OQ). 

Proposition 9. T/ie solution to (TT2J) determined by the A m ^ n through 
(20) 

/" 1 1 

ym,n — / n^m,^ yOy^ ) a - L to Sim,n C i/m,n 

It thus suffices to study ( TJgj) . 

Proof. Taking x = in ( T181) we obtain (1X9]) : using now ( {19]) in ( fLSi) we 

have 

1 .,-1 e -«m,n[*-* , [^ b ( a /)d ;E ' 



(21) y m>n — / ~^ K m, 



+ [2A min - (An+1,+1 + An+l,n-l)] 

— 1 p~ K m,n \x — x'\ I I I \ r jr r l _|_ p — Km,n \x\ A r ^ K m,n\x\ 

ci m,n </AH x it -"-1X1,11 e ym,n 

□ 

Remark 10. If y G then A m ^ n = y m>n (0) satisfies (TJJ). 
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Let = y^ n (0) where y^ nn is a solution to ( |T6l) . The solution 

of (|T6j) has the freedom of a multiplicative constant; we choose it by 
imposing 

(22) A° 00 = lim(<r - l)4,,o 

' CT— >1 

It is clear A° m n satisfies the homogeneous equation associated to ffT9l) 

(23) (Va + mj - imX - l)A° mn = -^A° m+l n+1 - ^A^ +l n _ 1 
3.4. Positions and residues of the poles. Define 

(24) (x = 1 - - u 

-UJ- 

To simplify notation we take u > 1 in which case o"o = 1. The general 
case is very similar. 
Denote 

(25) B := {ma; + mX + i : m G Z, n G Z, m < 0, |n| < |m|} 

Proposition 11. T/ie system [2&\) has nontrivial solutions in M\, iff 
a = ao(=l as discussed above). If a = 1, then the solution is a constant 
multiple of the vector A° m n given by 

{A° m n = m > and (m, n) + (0, 0) 
^m,n = m < and m < n < —m 
A° m>n = l (m,n) = (0,0) 

and obtained inductively from (23\) for all other (m,n). (Note that 
a = 1 is used crucially here since allows for the nonzero value of 

Proof. Let a = 1. By construction, A defined in Proposition [TT] satis- 
fies the recurrence and we only need to check Since 

4° 4- 4° 

j^Q ^■ L m+l,n+l ~ ■ rs -m+l,n— 1 

2(Vc" + ntu + imX — 1) 



and V a + ncu + imX — 1 7^ 0, we have 

2 m 



ALJ < C- 



rn .11 



proving the claim. 

Now, for any a, if there exists a nontrivial solution, then for some 
n , m we have A^ o mo 7^ 0. By f[2"3"j) . we have either 

(27) K_i, mo+ il > \ \Wa + n u + im Q X-l)\ ■ K o>m J 
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or 



(28) 



rao+l,m +l 



> - M a + n uj + im \ - 1)| • \A\ 



o 

n ,m 



It is easy to see that if in u + m A + i G B c or a 7^ 1, the above 
inequalities lead to 



for large m > (note that in these cases \J o + ncu + imA — 1 7^ 0), 
contradicting (J5]). 

Finally, if a = 1, then A is determined by Aq via the recurrence 
relation (1231) (note that A | B c = 0). This proves uniqueness (up to a 



Combining Proposition [8] and Proposition [TTJ we obtain the following 
result. 

Proposition 12. TTie solution ip(p) to equation (Tip zs analytic with 
respect to \/—ip, except for poles in B. 

Proof. Proposition [11] shows that (1231) has a solution A for cr g £>; by 
Proposition [S], A has singularities in B, and the conclusion follows from 
Proposition O □ 

So far we showed that the solution has possible singularities in B. 
To show that indeed if) has poles for generic initial conditions, we need 
the following result: 

Lemma 13. Let H be a Hilbert space. Let K(o~) : H —>■ H be compact, 
analytic in a and invertible in B(0, r) \{0} for some r > 0. Let vo(a) ^ 
Ran(I — K(0)) be analytic in o. If v(a) G H solves the equation (I — 
K(o~))v(a) = fo(cr), then v(a) is analytic in a in 5(0, r) \ {0} but 
singular at o = 0. 

Proof. By the Fredholm alternative, v (a) is analytic when o 7^ 0. If 
vio) is analytic at a = then t> is analytic and Vo(a) G Ran(I — K(0)) 
which is a contradiction. □ 

The operator €. is compact by Remark [TJ The inhomogeneity yo in 
equation (ITSl) is analytic in ^fa. Furthermore, at a = 1, Ran(I — (£) 
is of codimension 1 (Proposition [TTT) . Combining with Lemma [TBI we 
have 

Corollary 14. For a generic inhomogeneity y , y(cr) is singular at 
a = 1. Equivalently, ip{p) has a pole at p = i. 



(29) 




constant multiple) of the solution. 



□ 
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It can be shown that tp{p) has a pole at p = i for generic i/jq. We 
prefer to show the following result which has a shorter proof. 

Proposition 15. The residue i?o,o of the pole for 'if) at p = i is given 
by 

(30) #0,0 = lim(<x - l)A ,o = - A M + 2^ ,o] CT=1 

In particular, i?o,o 7^ /or /arge A and generic initial condition ipQ. 
Proof. When m = and n = ( Fl9i) gives 

(31) (v 7 ^ - 1)A 0)0 = -g^L-i - 2 Al ' 1 + 9o '° 

Clearly Aq^ is singular as a — > 1, which implies that ^ has a pole at 
p = i with residue given in fl30|) . Thus i?o,o is not zero if the quantity 
[-^4i,-i - A ±>1 +2g 0;0 ) a=1 is not zero. First, g ,oU=i is not zero by 
definition: 

S'o.ol^i = i / -e~^if> (a/)da/ 

./HP ^ 



Next, taking m — 1, n — 1, and cr = 1 in ( 1191) we obtain 

1 , 1 , 
---42,2 - -A 2 fi + 9lA 



(Vl +cj-«A- l)i4 M 



CT=1 



Thus for any c > when A is large enough we have 

IAi.il < cr 1 [|0 lfl | + max{|A 2)2 |, |A 2i0 |}] CT=1 

Estimating similarly v4 22 and A 2 ^ and so on, we see that \A lt i\ = 
0(c~ l ). When c is large enough we have \A lf i\ < |go,o|- Analogous 
bounds hold for Ai^i showing that \—Ax-i — + 2o ,o] CT= i * s no ^ 
zero. 

□ 

Corollary 16. For generic initial condition ip has simple poles in B, 
and their residues are given by R m , n = A^ nn . 

Proof. Take a small loop around o = and integrate equation ([TBI) 
along it. This gives a relation among R m ,n which is identical to (|23|) : 

(32) (Vcr + nu - im\ - l)R m ,n = ~Rm+l,n+l ~ ^+l,n-l 

Proposition [15] and (l22i) implies that i?o,o = ^oo- The rest of the proof 
is follows from Proposition [TTJ □ 
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Remark 17. It is easy to see that there exist initial conditions for 
which the solution has no poles. Indeed, if the solutions ipi an d ip2 
have a simple pole at p = i with residue a\ and 02 respectively, then for 
initial condition t/> ,o — Q^^o — a i4 ! 2,o ! the corresponding solution ip 
has no pole at p = i. 

3.5. Infinite sum representation of A m ^ n . Taking a = 1 in ({19]) we 

get 

1 . 1 



(33) (Vl + nuj - im\ - l)A m>n = --An+i,n-i - -A rn+ltn+1 + g m ^ n 

For r = (ai,...,ajv) G { — 1, 1}^, we define rj 1 = (a±, 0, 0). 
(Note that r = r£). We denote Erf = o< and {-1, 1}° = {0}. 
Let 

o 1 

-Dm 1 



\/T + nu — imX — 1 
and for some r G {—1, 1}^ define 



JV-l 



BrnriN — B mn ^{T) — ]TJ B m+ j n+ ^- T 

j=0 



Equation ( 1331) implies 



A m ,n — \~ -LJ / ., D mnN-l rt-m+Nji+Y^T 

(34) ^ x 

3=0 r6{-l,lp' 



As iV — > 00 we have 

N 



1 



j=0 



and A TO) „ goes to zero as m — > 00, and thus we have 

lim (— 1)^ — B mnN _i A m+Nn+ y' T — 

iV— >oo Z ^- ' 

re{-l,l} JV 

In the limit N — ► 00 we obtain 

00 j 

(35) A m ^ n ^ l^'^y ^ ^ B mn jg m j r j n j r ^ T 

j=o re{~l,l} 4 
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Remark 18. Truncating the infinite expansion to N, the error is 
bounded by 



(36) 



1 N 

TyN (II A™+i,n+£r° Am+N,n+J2r) 



2 N 

re{-l,lp 3=0 



4. Proof of Theorem Q] 

In §3.41 it was shown that for a generic initial condition ipo(x), the 
solution ^(xjp) has simple poles in B, with residues R m ,n = An,n- 

Since y G the inverse Laplace transform can be expressed us- 
ing Bromwich contour formula. Recall that y differs from the original 
vector form of ip by ( TTTT) . we have 



1 



(1) i/>(x,t) = C-^(x,p) = C-\h l ) + — / e pt M*,P)dp 



c+ioo 



c—too 



The fact that y G Jtf, also implies that ipi(x,p) — > fast enough as 
p — ► c ± zoo. Thus the contour of integration in the inverse Laplace 
transform can be pushed into the left half p-plane, after collecting the 
residues. As a result, for some small c < the contour becomes one 
coming from c — zoo, joining c — ie, 0, and c + ie (for arbitrarily small 
e > 0) in this order, then going towards c + zoo. 

Thus we have 



(2) 



il>(t,x) = + Res\ p=i {e pt ^ l ) 

+-^—e ct J e lst {^>i{x, c + is) +tp 1 (x,c- i s) j ds 
1 f c ~ ie „♦ ■> , , , 1 



c+ie 



e pt ijji(x,p)dp H / e pt *jj 1 (x,p)dp 

2m J 

By Corollary [16] we have 

Resl^e^) = #o,o = A° >0 

The third term in (jSJ) decays exponentially for large t (since the 
integral is bounded), while the last two terms yield an asymptotic power 
series in 1/y/i, as easily seen from Watson's Lemma. 

Combining these results and the fact C~ l (hi) = o(l/t) (Remark [6]) 
concludes the first part of Theorem [TJ with r(X,u) = Ro,o- The rest 
follows from Proposition! 
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5. Proof of Theorem [2] 
When lu — 0, the equation 

i ^ = ( - ^ - 25(x) + 26 (x) e- xt coa(ut)) tp 



dt \ dx 2 
becomes 

dip ( d 2 



26(x) + 26{x) e~ A ') ip 



dt V dx 2 
Rewriting A m n and g m>n as A n and g n , (fl9l) becomes 



(Vct - miA - 1) A n = -A n+ i + # n 
Since a; = 0, (155]) simplifies to 

OO I 

f± a/1 - in H 



(2 ) ^D-'i'n ,.,' 



j=o j=o a/I- '^^+J')^- 1' 

5.1. Proof of Theorem [2], (i). When n — 1 ([T]) becomes 

it 

With the notation 

A; 

^=n(v /i -^' A - 1 ) 

D becomes 

If 9k 
h k 

h k = e Wk ^J\ k - l {k-l)\ 

We have 

w k+1 -w k = log(Vl -ikX - 1) - - log(Afc) 



k=i 

Let 



Differentiating in A we obtain 



d , , I 



d\ y ^ 2AV1 - ikX 

Let Mfc be so that 

d d i \J — i k 

Tx Uk = Tx Wk ~ 
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Then, 



. . d . . i sj —i k — % i y—ik 1 

(3) Tx {Uk+l ~ Uk) = A^ + ^7^ + 2AVl-^A 

By taking the inverse Laplace transform of Q in k we get (we use p 
as the transformed variable here) 

Integrating (j4j) with respect to A gives 

VT\ (-2 + 2e- p - z 3 / 2 v^FAerf (^fi) 
jC. -Hi,- = 



2(-l + e-f)v^F(pA) 3 / 2 
Thus 

1 e~ Wfe 



/ifc v/A*- 1 ^- 1)! 

,oo _ v ^A(-2 + 2e-P-« 3 / 2 v ^¥Aerf(^^ 
= eXP (-/ c ~* 2(-l + e- P )0F(pA)3/ 2 ~* 

/ H\/—i k\ , j-fc 1 
x exp = — ) A 2 



Finally we obtain 



e -p /-c+ioo 22 V /Z ^\,i^ 1 

5f fc exp — — A 2 



poo 

exp 



v^A (-2 + 2e~* - z3/2 v ^lerf(=^) 
" P 2(-l + e- P )0F(pA)3/ 2 -*)^ 



5.2. Proof of Theorem (2J, (ii). Here we assume that expansion of 
A\ as A — > is invariant under a | rotation; that is, there are no 
Stokes lines in the fourth quadrant; this would be ensured by Borel 
summability of the expansion in A. 

Let A = ir with r < 0, and for simplicity let g = 1, then (|2j) implies 
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A" (-i)* f^ nLi(v / rT^+i) 

^ (-l) k k\r k 

n=l v ; 

The Euler-Maclaurin summation formula gives 

n rn 

Y^log(Vl + kr + 1) ~ / log( Vl + xr + l)dar + C 



(6) 



where 



fc=i 



1 , . , /- z— „, l , a/T+It 

— + A; log (VI + fcr + l - -k + - + C 

r 2 r 



l / r nl/r 

C ~ ^ log(Vl + kr + 1) - / log(vTTar + 
fc=i Jo 



ljaa; ~ 



Therefore 



^ exp(-l + fclog (yT+fcf + 1) - |fe + ^ ) 
(7) A x ~ 



(-l) k k\r k 
k=i K ' 



Since 



exp (_i + k bg(vr+g + 1) _ ik + 

(-l) fc A;!r fc 

exp ( _ A + + V> _ _ !2iM + 10 «(-'') 



2r 3 v r 2 2 2 

applying the Euler-Maclaurin summation formula again gives 

2 l/3 3 l/6 r( 2 )e -f ( _ a) l/6 



A 1 



2VvF 



5.3. Numerical results. Figure [T] shows log(|i?o|) as a function of 
log(A), very nearly a straight line with slope 1/6 (corresponding to the 
A 1 / 6 behavior), with good accuracy good even until A becomes as large 
as 1. 
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O 

o 



Figure 1. Log- log plot of \R \ function of A for 
L) — 0. Rofi is the residue of the pole of ip(p, x) at p — i, 
see Corollary [161 



6. Ionization rate under a short pulse 

We now consider a short pulse, with fixed total energy and fixed 
total number of oscillations. The corresponding Schrodinger equation 
is 

(1) z ^ = ( - - 26(x) + 2A 8(x) e~ xt cos (art)) ij> 

where A is now a large real parameter (note the A in front of the expo- 
nential). We are interested in the ionization rate as A — ► oo. 

By similar arguments as in §3.51 we have the convergent representa- 
tion 

oo f \ \ L * 

(2) A mjU = y^-iy i-j 2^ II Bn+j,m+\T9\9n+i,m+\T\ 

i=0 ^ ' T e2 i j=0 

Figs. [2] and [3] give |-Ro,o| (see Corollary [16]) and A under different u>/\ 
ratios. On small scales on the A axis, |-Ro,o| exhibits rapid oscillations, 
easier seen if uo/ A is smaller. 
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12 3 4 
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lambda 



Figure 2. |i?o,o| as a function of A, with fixed ratio 
u/\ = 20. i? ,o is the residue of the pole of ^}(p,x) at 
p = i, see Corollary [161 

7. Results for A = and to ^ 

We briefly go over the case A = 0, where ionization is complete; the 
full analysis is done in [TT| . In this case ip does not have poles on the 
imaginary line; we give a summary of the argument in [llj . 

The homogeneous equation now reads 

(1) \/a + mw A m = -~An+i - g^m-i + A m 



Thus we have 

^ \fa + mu A m A m = -- A rn+1 A m -^^2 A rn-\A m + ^ A r , 



The first sum and the second sum on the right hand side are conju- 
gate to each other, and each term in the third sum is real. So the right 
hand side is real, thus the left hand side is also real. 

For 3(a) 7^ 0, 3 (Ja + mu A m A m ) has same the sign as Scr. There- 
fore the sum can not be real and the equation has no nontrivial solution. 
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lambda 



Figure 3. |i?b,o| as a function of A when u/X = 5. _Ro,o 
is the residue of the pole of i/j(p, x) at p = i, see Corollary 
UM 

When 3(a) = 0, for m < 0, all 3 (a/ct + mu; A m A m ) have the same 
sign and for m > 0, a/<t + mw A m A m is real. Since the final sum is 
purely real, this means A m = for m < 0. But then, recursively, all 
A m should be 0. 

Zero is thus the only solution to ([jQ). By the Fredholm alternative 
the solution A is analytic in y/a and thus the associated y is analytic 
in yfo. This entails complete ionization. 

7.1. Small A behavior. We expect that the behavior of the system 
at A = is a limit of the one for small A. However, this limit is very 
singular, as the density of the poles in the left half plane goes to infinity 
as A — ► 0, only to become finite for A = 0. Nonetheless, given a A, small 
but not extremely small, formula (|35|) allows us to calculate the residue. 

Figure H] shows the behavior of the residue versus u, for A = 0.01. 

We show for comparison the corresponding result when A = 0, in 
Figure [5j 
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Figure 4. |-Ro,o|, at A = 0.01, as a function of uj. i?o,o 
is the residue of the pole of x) at p — i, see Corollary 




r=0.01 



CO 

1/51/4 1/3 1/2 1 

Figure 5. log 10 r _1 , at A = 0, as a function of uj/ojq. 
T is the Fermi Golden Rule exponent for the probability 
decay and hu Q = E , the energy of the bound state of 
one delta function of amplitude r [7]. 
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